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ABSTRACT 

We have obtained Keck HIRES spectra of three images of the quadruply 
gravitationally lensed QSO 2237+0305 to study low ionization absorption systems 
and their differences in terms of projected velocity and column density across 
the lines of sight. We detect Call absorption from our Galaxy, and a system 
of High Velocity Clouds from the lensing galaxy (z=0.039) with multiple Call 
components in all three sightlines. Unlike the situation in our Galaxy there is no 
prominent Call absorption component (with an equivalent width exceeding 60-70 
mA) close to the velocity centroid of the lensing galaxy 2237+0305. Instead, Call 
components with total equivalent widths similar to those of Galactic intermediate 
and high velocity clouds are spread out over several hundred kms~^ in projection 
along the sightlines at impact parameters of less than one kpc through the bulge 
of the galaxy. A Call absorbing thick disk like in our Galaxy does not seem 
to extend into the bulge region of the 2237+0305 galaxy, whereas high velocity 
clouds seem to be a more universal feature. We have also studied three low 
ionization Mgll-Fell systems in detail. All three Mgll systems cover all three 
lines of sight, suggesting that the gaseous structures giving rise to Mgll complexes 
are larger than ~ 0.5 kpc. However, in most cases it is difficult to trace individual 
Mgll 'cloudlets' over distances larger than 200-300 h^Q pc, indicating that typical 
sizes of the Mgll cloudlets are smaller than the sizes inferred earlier for the 
individual clouds of high ionization gas seen in CIV absorption. We tentatively 
interpret the absorption pattern of the strongest Mgll system in terms of an 
expanding bubble or galactic wind and show that the possible loci occupied by the 
model bubble in radius-velocity space overlap with the observed characteristics 
of Galactic supershells. 
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Subject headings: galaxiesrISM — galaxies: bulges — ISM:bubbles — quasars: 
absorption lines - gravitational lensing - quasars: individual (Q2237+0305) 

1. Introduction 

Absorption lines in the spectra of the separate images of gravitationally lensed QSOs 
provide a unique opportunity to investigate the detailed structure of the interstellar gas and 
the intergalactic medium on scales of ten parsecs to a few tens of kpc. The ultimate goal of 
such observations is to study the evohition of the microstructure of the gas in galaxies and 
in the intergalactic medium from the epoch of formation of the first luminous objects in the 
universe down to the present epoch. For the past few years we have undertaken a survey 
of lensed QSOs with the HIRES spectrograph (Vogt et al 1994) on the Keck I telescope on 
Mauna Kea. The first paper in the present series (Ranch, Sargent and Barlow 1999) analyzed 
a low ionization absorption system at Zabs = 3.538 on a scale of only 26h^Q^ parsecs. The 
second paper (Ranch, Sargent and Barlow 2001) investigated the structure of C IV absorbing 
clouds, including an estimate of the rate of input of mechanical energy, while the third paper 
(Ranch, Sargent, Barlow and Carswell 2001) set a hmit on inhomogeneities in the "Lyman 
alpha forest" clouds at redshift z ~ 3.3. 

In the present paper we turn our attention to low ionization systems which are charac- 
terized by Mgll or Fell absorption lines. This gas is likely to be denser, more metal-rich, 
and more closely associated with galaxies than either the average Lyman a forest or the high 
ionization CIV absorption systems. 

At intermediate redshifts Mgll absorbers have usually been interpreted as the large 
gaseous halos of luminous galaxies. The sizes of entire Mgll complexes are relatively well 

constrained: Smette et al (1995), from observations of HE1104-1805A,B, estimate lower 
limits of > 20h^Q kpc for the diameters of Mgll systems (see also the high resolution study 
of the z=1.66 system in that QSO by Lopez et al 1999). Monier et al (1998) find an upper 
limit of 20/i^Q^ kpc for a strong Lyman limit absorber. HIRES spectra we have obtained of 
the two lines of sight to UM673 A and B (to be published) show three strong Mgll absorption 
systems. The one at z=0.426 (transverse separation between the lines of sight: lA.dh^Q kpc) 
shows absorption in only one of the images, whereas in the other two systems (z=0.492, 
sep. = 15.5h^Q kpc; z=0.564, sep. = 13.9/i^q^ kpc) there is Mgll absorption of similar 
strength in both sight lines. Siirveys for galaxies in the vicinity of Mgll absorption systems 
and attempts to reconcile the total absorption cross section with the luminosity function of 
galaxies (Bergeron & Boisse 1991; Steidel 1995; Bowen, Blades & Pettini 1995; Churchill et 
al 1999) have shown that galaxies appear to be surrounded by halos of Mgll gas with radii 
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estimated to be on the order of 60-130 h^^ kpc. 

The spatial extent of individual Mgll absorption components as opposed to entire Mgll 
systems is somewhat harder to pin down. The low ionization system discussed by Ranch, 
Sargent and Barlow (1999) showed that Sill and CII gas cloudlets (at least at z=3.6) can 
be as small as a few tens of parsecs. Petitjean et al (2000) have investigated intermedi- 
ate redshift Mgll systems in a (spatially unresolved) HIRES spectrum of the lensed BAL 
APM08279+5255 and conclude from the residual intensity in the Mgll doublets that the 
sizes of individual clouds are less than I.S/ijq^ kpc. The only other independent size esti- 
mate on individual cloudlets comes from ionization calculations (e.g., Bergeron & Stasinska 
1986; Steidel 1990), which, however, provide more indirect constraints. Rigby et al (2002) 
using such arguments, find that a subset of iron-rich Mgll systems may be consistent with 
cloudlets as small as 10 pc in diameter. 

In the present paper we address the question of the small scale structure of low and 
intermediate ionization gas in more detail. In particular, we investigate the absorption lines 
in the spectrum of the "Einstein Cross" in which Q2237-I-0305 (zgm — 1-69) is gravitationally 
lensed by a 15th magnitude galaxy with a redshift z = 0.0390 (Huchra et al. 1985). Subse- 
quent work (Yee 1988) revealed that the galaxy is a barred Sb with a ring, and has a disk 
scale length tq = 9.4" and an effective radius for the bulge Ve — 1.9". Q2237+0305 turned 
out to have four images roughly arranged in the form of a cross centered on the nucleus of 
the galaxy (Yee 1988; see his fig. 4). The typical separation between images is less than 2". 
Crane et al. (1991) obtained B = 17.96 magn. for component A, B = 17.82 magn. for B, 
18.66 magn. for C and 18.98 magn. for image D. Later observations showed that the relative 
brightnesses of the images is subject to changes, in part due to microlensing (Irwin et al 
1989). Note that the flux ratios obtained during the current observation (in the caption of 
Fig. 1) are very different indeed. 

The original discovery spectra only showed self absorption in the Lyman alpha and C IV 
emission lines. Hintzen et al. (1990) searched for but failed to flnd Ca II H and K absorption 
due to the lensing galaxy in the spectra of the quasar images. However, they detected Mg 
II and Fe II absorption at Zabs — 0.5664 and a C IV complex at Zahs — 1-694 and 1.697. 

In section 2 we describe the high resolution observations of the three brighter components 
of Q2237-I-0305 with HIRES. Section 3 presents the detection of Call absorption from the 
Milky Way and Section 4 describes a Call system arising in the lensing galaxy. Section 5 
discusses the properties of a pecuhar, strong Mgll absorption and tries to give a tentative 
interpretation. Sections 6 and 7 briefly describe the other two low ionization systems towards 
Q2237-I-0305, and section 8 gives the conclusions. 
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2. Observations and data analysis 

Q2237+0305 has four images roughly arranged in the form of a cross centered on the 
nucleus of a 15th magnitude galaxy. A picture of the lensed QSO taken off the Keck HIRES 
guider is shown in Fig. 1 where the four images are identified. Spectra of the three brighter 
images A, B and C were obtained with HIRES in October 1998. The blue cross disperser was 
used, giving a wavelength coverage of 3645-5204 A which puts the C HI] 1909 emission hne 
at the red end of the range covered. We used a 0.86" slit, giving a resolution of 6.6 km s"^. 
The exposures, broken into 3000 second sections, totaled 18,000 seconds on image A, 18,000 
s on B and 15,000 s on C. A heliocentric correction was applied. All wavelengths quoted here 
are vacuum values. The individual spectra of each image were combined and scaled to the 
same continuum level in order to compare the details of the absorption lines. The absorption 
systems were fitted with the Voigt profile fitting routine VPFIT (Carswell et al 1991, 1992). 
Fitting components were put in until a statistically satisfying value was achieved. Usually 
this involved obtaining a reduced with a probability > 1% to have ocurred by chance. 
The continuum level was treated as an additional free parameter (see paper HI) to reduce 
systematic errors that could have crept in with the interactive continuum fitting. 

The interpretation of this dataset is somewhat more challenging than usual because the 
QSO images closely surround (to within an arcsecond) the peaked emission from the bulge 
of the foreground lensing galaxy, so there is a high background level which rapidly changes 
in the radial direction (Fig. 1). The data show several instances of absorption lines which to 
the eye appear fiat-bottomed and saturated while not dropping down to zero intensity in the 
line center. This phenomenon is most obvious in the strong absorption of the Mgll complex 
at z=0.5656, especially in images B and C (see Figs. 5 and 6 and below, and the discussion 
in section 4.1), where there appear to be partly flat-bottomed Mgll 2796 lines with a residual 
intensity of about 11% (in image B). At values of 1.35 and 1.16 the Mgll doublet ratios for 
the strongest components in the A and C images, respectively, also suggest partial saturation. 
There are several possible explanations for this behaviour: (a) the absorbing gas cloud does 
not cover the continuum source of the background QSO and the troughs are partly filled in 
with light by-passing the absorber; (b) the absorption lines are ver?/ narrow and unresolved by 
the instrument; and (c) there is a foreground light source which partly fills in the absorption 
lines of systems arising between the foreground source and the background QSO. Possibility 
(a) is quite unlikely in that all of the low ionization systems are very far away from the QSO 
redshift. Case (b), unresolved narrow lines, may indeed by present, but the fiat bottoms of 
the absorption lines and the fact that different lines often reach similar levels argues against 
it being the main effect: it would require a conspiracy between individual optical depths 
and velocity separation between the components to produce such a pattern. Possibility (c), 
pollution by the light of the bright lensing galaxy appears most likely: inspection of various 
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strong lines in all three lines of sight shows that the residual intensities at the line bottoms 
correlate with the local flux of the QSO (they are weaker where the QSO is brighter, e.g., in 
the CIV emission line), as one would expect if the contribution of the QSO continuum light 
to the total spectrum increases relative to the light of the foreground galaxy. The relative 
contribution of galaxy and QSO continuum to the flux at a given wavelength, Fgai/ FqgQ, 
can be written in terms of the residual intensity at the bottom of an apparently saturated 
absorption line, F^^*(A), as 



For image A the Mgll 2796 line at z=0.5656 gives approximately Fgai/ Fq^q ~ 0.08, while the 
strongest CIV 1548 line at z=1.693 gives 0.04. The CIV z=1.693 line is an associated system 
riding almost on top of the CIV emission line and hence may be a candidate for partial 
continuum coverage. However, we ignore this effect here since it can only make the line less 
black. For image B the Mgll 2796 at z=0.5656 gives approximately Fgai/FQgQ ~ 0.11, while 
the CIV 1548 hne at z=1.693 gives again about 0.04. In other words the contribution of the 
hght from the lensing galaxy to the integrated spectrum is up to about 11% near the Mgll 
2796 z=0.5656 system. 

A naive measurement of the column densities without correcting for the light pollution 
by the lens would underestimate the column density for weak lines on the linear part of 
the curve of growth by on the order of - 12% (see section 4 below), but for saturated 
lines they may be highly underestimated. Therefore, the column densities and measurement 
uncertainties in tables 1,2, 4 and 5 below (which reflect only the proflle flts without correction 
for the residual flux from the lens) should be taken as lower limits. The exception is the 
z = 0.5656 low ionization system (section 4) where many lines are obviously saturated, so 
a correction was attempted by simultaneously adjusting the zero level of the spectroscopic 
flux while fltting the absorption line proflles. The correction was possible because there 
were enough saturated regions to constrain the zero flux level, but it could not be applied 
meaningfully to the other absorption systems which appear to be mostly optically thin and 
do not have enough transitions available. Thus table 3 contains actual estimates of the 
absorption line properties. 




(A) l-F'-^'iX) 



F'-^'iX) 
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3. Call absorption from the Milky Way 



The Call absorption from our Galaxy can be modelled satisfactorily with a single veloc- 
ity component in all three images (see the flt parameters in table 1). A pairwise comparison 
between any two of the three sightlines is shown in Fig. 2. The column densities and Doppler 
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parameters are identical to within the errors, and the projected velocities along the line of 
sight differ by less than 2.9 a between the images. In other words, there is no significant 
difference in the appearance of the Galactic Call systems in the three lines of sight. The 
mean weighted heliocentric velocity of the three sightlines is cz = —13.9 ± 0.5 kms^^. The 
direction of the Q2237+0305 sight lines in galactic coordinates (/^^ = 71.8, b^^ = -46.1 deg) 
subtends an angle of ~ 63 deg with the heliocentric velocity vector of the local bulk flow 
of interstellar cloudlets within ~ 30 pc (wheZio = —28.1 ±4.5 kms~^; Prisch, Grodnicki, & 
Welty 2002) so the projected velocity of that flow relative to the 2237+0305 line of sight 
would be -12.9 kms~^. This is well consistent with the idea that the Call observed here 
arises from the cluster of local interstellar gas clouds within 30 pc. The maximum angle be- 
tween the three lines of sight to Q2237+0305 is about 1.8 arcsec. Thus, if the Galactic Call 
absorber towards 22327+0305 resides within 30pc, our observation means that there are no 
significant detectable differences in that particular Call cloudlet over transverse separations 
smaller than 54 astronomical units. 



4. High velocity CeJI clouds related to the lensing galsixy (z=0.039) 

4.1. Observed properties 

We detect an absorption system at z=0.039, the redshift of the lensing galaxy 2237+0305. 
The system is visible in all three fines of sight and consists of moderately narrow (7.6 < b < 
24.3 kms~^), presumably nonstellar Ca II H and K absorption fines. A comparison of the 
Ca II H and K profiles in each pair of the three images is shown in Fig. 3. Each line of sight 
is shown separately in Fig.4.^ Any underlying, broad stellar Ca II absorption feature due to 
the lensing galaxy has been removed during the background subtraction and by drawing the 
continuum such as to eliminate large scale features. The velocity scale in Figs. 3 and 4 is 
relative to the lensing galaxy redshift Zabs = 11696(±40) kms~^/c = 0.039014 given by Foltz 
et al (1992). 

Fig. 3 shows that image A has a prominent absorption feature at +103 km and a 
weaker feature at -91 km s^^. ^ Image B has three (and possibly more) weak components 
at -288, -202, and +105 kms""*^. Image C has two strong components at -153 and -95 kms~"^ 
and two additional weaker ones at +38 and +90 kms^^. Only the two strongest lines, the 



"^Note that the wavelengths of Ca II H and K quoted in Fig. 2,3 and 4 are the vacuum values. 

■''To avoid confusion we do not show the fits on top of the spectra but give the fitting parameters redshift, 
Doppler parameter h and logarithmic column density log N and their 1 — a errors in table 2. 
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+103 kms~^ component in image A and the -95 kms~^ component in C are strong enough 
to be clearly present in both transitions of the Call doublet, and it is only for those that we 
can be sure of them being Call. Any of the other lines could in principle individually be an 
unidentified interloper from other transitions (there is one such case, marked by an x in the 
Call 3970 regions in Fig. 3 which is the 1550A line of a stray CIV system at z= 1.661), but 
the line density is much too low for this to explain the majority of components. 

Using intermediate resolution spectroscopy of the integrated light from the four QSO 
images, Hintzen et al (1990) have searched for Call from the lensing galaxy, deriving a 3 — a 
upper limit of 72 mA on the rest frame equivalent width of the Call K line. They compared 
this limit with the equivalent width of Call from the thick disk in our Galaxy and concluded 
that the absence of strong Call rules out the presence of a thick disk near the center of 
2237+0305. Calculating the Call 3935 A equivalent widths from the column densities given 
in table 2 we get W = 97+4, 122+12, and 424+9 ml, for the total equivalent widths of 
all the fitted components in images A, B, and C, respectively. The unweighted mean of the 
three images is 214 mA. These values are not corrected for the zero level problem mentioned 
above but if anything the real values can only be larger (probably by not more than 12%). 
The equivalent widths are formally inconsistent with the null detection of Hintzen et al. 
However, their result is an average of the Call absorption along all four lines of sight which 
they observed through a 2.5" diameter circular aperture, whereas our observation was done 
only for the brightest three images individually with a slit 0.86" wide. The fourth line of 
sight (included in their average but not in ours) may not have any absorption, which would 
reduce the discrepancy. Another possible explanation has to do with the uncertainties due 
to pollution of the QSO spectra by the bright foreground galaxy. Hintzen et al used an 
aperture larger by more than a factor 8 than ours (assuming that we had 0.7" seeing), and 
the depth of the stellar H and K lines in their spectrum shows that the galactic contribution 
to the QSO spectrum was severe. Hintzen et al dealt with this problem by subtracting 
sky "from a position 36 arcsec offset from the QSO", and "subtracting a suitably scaled 
spectrum of the M31 globular cluster K225-B280" to remove the broad stellar H and K 
lines. In our case, the sky background (including the stellar contribution) was subtracted 
in windows immediately adjacent to the QSO image. We can only speculate that given our 
own problems in subtracting the spatially fiuctuating sky background (discussed in section 
2) it is conceivable that Hintzen et al may have had similar ones. They may have subtracted 
less of the stellar background than we did which would have led to them underestimating 
the depth of the absorption and thus the Call equivalent width. 

In any case. Call is detected and it appears that the measured total equivalent width 
is not dramatically different from what is seen in the Galaxy: Bowen (1991) found the 
typical high-latitude Call equivalent width distribution to peak at Wsinb = 110ml for 
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lines of sight going outward from the Galactic disk. In the present case the lensing galaxy 
at z=0.039 is inclined by 60 degrees (Irwin et al 1989) with respect to the plane of the 
sky so the equivalent width expected based on the Bowen study would he W = 220mA . 
Our measurement gives the two-sided equivalent width (intersecting the whole galaxy) and 
translates into 214/2 = 107 mA for a line of sight going out from the central plane of the 
galaxy. This is about half of what is seen in the Milky Way for the thick disk. However, most 
of the absorption components in the HIRES data formally qualify as intermediate or high 
velocity clouds (e.g. Wakker 2001), so a comparison of their equivalent width with a subset 
of higher velocity clouds in our Galaxy would be more appropriate. Galactic intermediate 
velocity clouds studied by Morton & Blades (1986) lead to an average of 100 ± 50 mA 
(Hintzen et al. 1990), similar to the present result, 107 mA . 

We do not have information on the HI column densities so we cannot obtain the gas 
phase abundances for these clouds. If the abundances measured in the interstellar medium 
of the Galaxy were similar to those in the 2237+0305 galaxy, then one can use the Call- 
HI abundance relation of Wakker & Mathis (2000) to determine the HI column densities. 
According to their formula, the Call column densities in table 2 would roughly correspond to 
logarithmic HI column densities between 18.6 and 20.1, with the exception of the strongest 
component in the C image (at -95 kms^^). This one formally gives a very high HI column 
(23.1) but it is beyond the range of validity of the Wakker-Mathis relation, and may be more 
indicative of a higher than usual Call gas phase abundance. 



4.2. Origin of the Call velocity structure 

The velocity structure of the Call system is different from that expected from absorption 
arising in a thick disk: 

The center of Call absorption for all Call components in all three sightlines (defined 
as the column density weighted mean velocity relative to the redshift of the lensing galaxy) 
occurs at -71 kms~^ with respect to the Foltz et al (1992) velocity of recession, 11696(±40) 
kms^^, or at +4 kms^^ relative to the HI 21cm velocity, 11621(±40) kms"^ given by Barnes 
et al (1999). Thus, within a 2-sigma deviation the column density weighted velocity centroid 
of the Call absorption agrees with the centroid of the stellar and HI mass distribution. 
However, Call absorption from the disk of our Galaxy seen against high latitude stars tends 
to have a strong component within about 10-20 kms~^ of the local standard of rest (e.g., 
Greenstein 1968; Blades & Morton 1983; Songaila et al 1986; D'Odorico et al 1989; Robertson 
et al 1991; Meyer & Roth 1991; Ho & Filippenko 1995; this paper). Similar Call absorption 
near systemic galactic velocities is also seen towards other galaxies, e.g., against SN1989M, 
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SN1993J, and SN1994I in the disks of M81 (Vladilo et al. 1994), NGC4579 (Steidel, Rich 
& McCarthy 1990), and NGC5194 (Ho & Fihppenko 1995), respectively, and towards the 
nucleus of M31 (Morton & Andereck 1976). This is in contrast to the present case where 
there does not seem to be a zero velocity component in any of the three sightlines, no matter 
which of the above velocity centroids for the galaxy is adopted. 

Furthermore, the spread in velocity of the Call components is larger than expected if 
caused by the disk. The total unweighted velocity widths (highest minus lowest velocity 
component) are Af^ = 193, Af^ = 393, and Afc = 243 kms^-*^ for the three images A, B, 
and C, respectively. Even if only the relatively certain strong components at -153 kms~^ (C 
image) and +103 kms^^ (A image) are used the velocity difference is at least 256 kms^^. 
The typical line of sight from the source QSO has an impact parameter of ~ 1" or l.lhgQ^ 
kpc. Therefore, the three lines of sight pass through the bulge of the lensing galaxy, rather 
than the disk, and it is clear that the rotation of the disk cannot be responsible for the 
large velocity widths even if it extended to the very center of the galaxy. Thus we agree 
with Hintzen et al, albeit for different reasons, that this case does not look like thick disk 
absorption. 

The origin of the large velocity width of the Call absorption is not obvious. The width of 

the absorption system is not unlike what is seen in sightlines towards the Magellanic clouds or 
in the Magellanic stream (Songaila 1981; West et al. 1985; Songaila et al 1986). It is possible 
that a similar configuration involving tidal gas fiows (sec also Bowen et al. 1994) exists here 
as well. In their search for HI emission in the vicinity of the lensing galaxy, Barnes et al 
(1999) have detected a group of galaxies which contains the lens. The two objects nearest 
to the lens (in projection), the dwarf irregulars NWl and NW2 (in their nomenclature), are 
separated from it by about 150 kpc and 125 kpc in the plane of the sky, and in velocity 
of recession by -380 and -264 kms~^, respectively. It is conceivable that this group may 
produce a tidal feature like the Magellanic stream, but the negative velocity differences of 
the closest objects with respect to the lensing galaxy 2237+0305 may be inconsistent with 
the positive velocities of some of the Call components. More suitable but fainter satellites 
of the 2237+0305 galaxy may have eluded detection. 

Even with the lack of evidence for a tidal mechanism there still remains a wealth of 
possible causes as have been discussed in the general context of Galactic high velocity clouds 
(which we do not wish to enumerate here; for a review sec Wakker, van Wocrden & Gibson 
1999, and references therein). One possibility not yet mentioned is that the Call clouds may 
arise directly in the bulge of the z=0.039 galaxy. Barnes et al (1999) have argued that the 
stellar velocity dispersion measured by Foltz et al (1992), corrected for observational effects, 
agrees with various mass models in giving a line of sight velocity dispersion ^ 145 kms~^ for 
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the bulge. The velocity differences between Call components quoted above exceed this value 
by at least 110 and possibly as much as 250 kms~^ (if all Call components arc real). To 
account for the additional velocity dispersion one could cnvokc outflows (winds, SN bubbles) 
from the bulge, or there could be other hydrodynamic effects in the bulge that stir up the 
gas of an otherwise quiescent disk. Morton & Andereck (1976) have observed additional 
velocity components at -260 and -450 kms^^ against the nucleus of M31 and argue that they 
are likely to be caused by outflows. However, there is no evidence for nuclear activity in the 
galaxy 2237+0305. 

Finally there are the size constraints from the three lines of sight. While the column 
density of the components is rather different between the sightlines there is a hint that 
at least two components of the C image (-95 and -|-90 kms~^ seem to correspond to weak 
features in the A image (at -91 and -1-103 kms~^). Accordingly, this particular system does 
have at least some coherence over 1.4 kpc, whereas the average cloud size (if defined as the 
distance over which the column density varies by 50% ) may be much smaller. Coherence 
on kpc scales may argue for clouds outside the bulge or disk, but the evidence is weak in the 
present case. Summarizing, the appearance of the z=0.039 system is consistent with that of 
high velocity clouds in the Galaxy. The velocity structure of the Call components and the 
apparent hole in any Call absorption by a thick disk argues for a kinematic connection of 
the Call cloudlets with the bulge and / or locations above the plane of the lensing galaxy. 



5. The Mgll-Fell system at z — 0.5656: expanding shell, collapsing sheet ? 

5.1. Observational properties 

This system is the strongest absorption complex in the spectral region covered by our 
observations (Figs. 5 and 6). The appearance of the Mgll 2796, Mgl 2853, and Fell 2600 
absorption lines makes it likely that the structure is a strong Lyman limit or damped Lya sys- 
tem with the hydrogen predominantly neutral: first, Mgl is rather strong (e.g., log A^(Mgl) - 
log A^(Mgll) = -1.42d=0.08 for the reddest Mgl/Mgll component in the A image (z=0.5669)). 
Photoionization models with CLOUDY (assuming an a = —1 powerlaw ionizing background 
with intensity J(912) = 10~^^ at z = 2.5, blueshifted to z=0.56) show that for a strong 
Lyman limit system (N(H1) = 10^^ cm~^) the total gas density would have to be rather high 
(> 350 cm~^) to get a Mgl/Mgll ratio as high as observed. For a damped Lyman a system 
with N(HI)=10^° cm~^, the density would need to be only as high as 0.5 cm~^. Second, 
strong Fell and Mgll hues often indicate damped systems (Boisse et al 1998; Rao & Turn- 
shek 2000) .The total equivalent widths of Mgll 2796 and Fell 2600 for the z=0.5656 system 
in hne of sight A are 1.20 A and 0.73 A respectively (after correction for the zero level). 
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According to Rao & Turnshek, 50% of those systems with EW(MgII 2796) > 0.5A and with 
EW(FeII 2600) > 0.5 A are damped (N(HI)> 2 x 10^° cm'^). 

The z — 0.5656 absorption system shows an interesting degree of coherence between 
the three hnes of sight (Figs. 5,6): in Mgll 2796, two massive groups of components (near 
130 and 175 kms~^ in the A spectrum (bottom)), separated by a gap in absorption near 
150 kms"'^, can be traced in all three images (with the A spectrum having an additional 
strong component near 255 kms^^). While the gap between the components does not shift 
by more than a few kms^^ between the spectra the two components seem to get further 
apart in velocity space as one goes from image A to C to B, i.e., in the order of increasing 
distance from A. The velocity differences va, vc, and vb, between the components are 
va — 69.8 kms"^, vc = 81.4 kms~^, and vb — 102.5 kms~^). The two massive components 
are clearly composed of multiple components indicative of differential motion in the gas; in 
particular the B spectrum shows a pronounced ragged appearance. The projected relative 
positions of the images A, B, and C, can be expressed in terms of distances AB (=0.66 kpc) 
and BC (=0.50 kpc) between images A and S,and B and C, respectively, and the angle 
Z{ABC){^ 47.5 deg). 

The fits to the absorption line profiles are hsted in table 3. Experience shows that Mgll, 
Fell and sometimes Mgl have a very similar component structure in velocity space even if 
the relative column density levels differ. Thus for each fitting component the redshifts of 
Mgll and Fell (and in the case of image A also Mgl) were tied together, i.e. they were forced 
to have the same redshift. The rcdshift errors arc quoted in the entry for the Mgll ion in 
table 3. Where the corresponding other ions at the same redshift have nominal errors of zero 
it should be understood that the Mgll redshift error is actually the combined error of the 
tied ions. 

To correct for residual fiux from the lensing galaxy the zero level within the fitting regions 
was adjusted as an independent variable simultaneously with the profile fits. The required 
adjustments of the zero level relative to a unit continuum in the three regions containing 
the Mgll doublet and Mgl (region 1), the Fell 2383 fine (region 2), and the Fell 2600 fine 
(region 3), turned out to be 7.9 ± 1.2%, -0.02 ± 4.8%, and -0.07 ± 7.9%, respectively, in 
the A spectrum, 10.6 ± 1.4%, 12.2 ± 4.8%, and 5.5 ± 6.6% in the B image, and 11.6 ± 0.7%, 
9.5 ± 1.4% and 10.7 ± 1.4% in the C image. These values can be used as F^^'^{X) in eqn. 1 
to obtain the ratio between galaxy and QSO fluxes which made it through the spectrograph 
slit. 

Limits on the gas temperature can be estimated from the line widths of the Mgll doublet 
and the Fell 2600 A line. An upper limit on the thermal width is given by the total line 
widths of the Mgll lines in Fig. 6. Limiting ourselves to the B image (which appears to have 
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the weakest lines, least likely to be saturated) a simultaneous fit to the Mgll doublet gives 
a weighted mean Doppler parameter 5.3 ± 0.4 kms~^ for the components in the z = 0.5656 
complex, corresponding to an upper limit on the temperature T < 4.1 x lO'* K. A more 
stringent limit may be obtained from comparison between two ions with different atomic 
masses, Mgll and Fell. For example, in Figs. 5 and 6 the Mgll 2796 spectrum of the 
B image shows a group of 3 - 4 ragged components between 140 and 210 kms~^. The 
three strongest Mgll components are well matched by corresponding Fell 2600 A lines (see 
also Fig. 5). Another good case for an unblended hne occurs near 250 kms~^ in the A 
image. Using those 3 Mgll/Fell components from image B where 6(MgII)> 6(FeII) together 
with the abovementioned line from the A spectrum, we employ the method discussed by 
Ranch et al (1996) to disentangle thermal and turbulent contributions to the total Doppler 
width. We have omitted two other cases where Fell was broader than Mgll as not providing 
information. The mean weighted thermal Mgll Doppler parameter turns out to be 4.0 ±0.8 
kms~^, corresponding to a temperature 2.3 x 10'^ K. Thus there appears to be a substantial 
thermal contribution to the broadening (on the order of a a few times 10^ K as opposed 
to a few hundred Kelvin or less which would not show up observationally) . Unfortunately, 
the large errors and the line widths being close to the resolution limit preclude any further 
constraints. 



5.2. Modelling the absorbing structure as a moving shell of gas 

It is impossible to say with certainty how this particular absorption system arises. How- 
ever, the double component structure mentioned above, the increase of velocity separation 
between the two components with transverse separation between the lines of sight, and the 
large coherence length (Fig. 6) are suggestive of the three hnes of sight intersecting two 
walls of a coherent, expanding (or collapsing) "shell" (Fig. 12). Earlier we described a small 
low ionization absorption structure in a high redshift galaxy towards Q1422+231 which had 
several of the absorption characteristics expected from an expanding shell (paper 1). The 
size of that structure, however, was at least one to two orders of magnitude smaller than the 
current object, which must be larger than several hundred parsecs in order to cover all three 
hnes of sight. It is conceivable that this absorption system could be associated with galactic 
HI shells and supershells which are known to be at least that large (e.g., Heiles 1976,1979; 
Deul & den Hartog 1990; Kennicutt et al 1994; Oey & Clarke 1997). Such structures may 
be caused by SN explosions or wind blown bubbles (Heiles 1979) or by the infall of gas from 
the intergalactic medium e.g., in the form of high velocity clouds (Oort 1967). Recently, 
Bond et al (2001a,b), in a detailed analysis, have given plausible arguments for at least some 
Mgll absorbers being caused by galactic ejecta, i.e., super-bubbles and -winds. The line of 
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sight component structure of the z — 0.5656 system is quite similar to that of several Mgll 
systems in Bond et al. 2001b. 

Here we attempt to model the absorption system as a spherical bubble in isotropic 
expansion (or contraction) . Having the velocity information from three lines of sight we can 
infer the radius and expansion velocity of the bubble from observations of the same system 
intersected by multiple lines of sight. The measured input parameters for such a calculation 
are the three velocity differences between the main components and the three transverse 
distances between the lines of sight. Unfortunately, these constraints are not enough to 
determine both radius R and expansion velocity independently, and only the relation ) 
can be obtained. The line of sight to the fourth [D) image also would have to be observed, 
which was not attempted by us because of the small separation and the faintness of image 
D. The Appendix shows how to calculate the relation R{v\ subject to the observational 
constraints. The result is given in the form of a diagram in Fig. 11, and the geometry is 
shown in Fig. 12. There is a broad minimum in the expansion velocity leading to values 
below 100 kms~^ for bubble radii between a few hundred parsecs and a few kiloparsecs, with 
the minimum expansion velocity, t; = 55 kms^^, occuring at a radius R = 0.83 kpc. At 
very small radii the diameter of the bubble approaches the transverse separation between 
the lines of sight. Such small bubbles require not only an increasingly unlikely alignment 
between the bubble walls and the lines of sight to still cover all three of them. They also 
lead to a near-grazing incidence of the lines of sight with respect to the bubble and they 
require very large expansion velocities to reproduce the velocity splitting seen between the 
components. At the opposite end, for very large bubbles with R » 1 kpc, the expansion 
velocity also would have to be very high: the spatial curvature of the bubble wall between the 
lines of sight would be smaller with increasing radius, and the differential velocity splitting 
between the major components seen in Fig. 6 can only be reproduced with an ever increasing 
expansion velocity. 

Galactic and local extragalactic superbubbles tend to show maximum sizes of about 1-2 
kpcs, and expansion velocities are typically below 150 kms~^ (Heiles 1979; Tenorio-Tagle 
& Bodenheimer 1988; Martin 1998). Thus the parameter space of observed superbubbles 
shows considerable overlap with the allowed combinations of R and v for our model (hatched 
area in Fig. 11). The observed total Mgll column densities in the two main components are 
2.7 X 10^^ and 3.3 x 10^^ cm""^ for the blue and red component in the A image, 1.8 x 10^^ 
and 2.7 x 10^^ cm~^ for B, and 7.8 x 10^^ and 4.5 x 10^^ cm~^ for the C image, respectively. 
We can obtain a crude upper limit on the size of the bubble by assuming that the absorption 
is arising entirely in a bubble shell composed of gas swept up by a wind outside the disk's 
interstellar medium (assuming subsequent cooling to produce Mgll as the dominant Mg 
ionization state). Then to produce a Mgll column density as strong as observed the gas 
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would have to come from a volume with radius 

R^3— ^2kvc( \ f ^tot ^Mg Z/Zq fMgIl \ ^ 

n ^ Vl0^''cm-2y VlO-2cm-3 3.3x10-5 0.3 0.5 / ' 

Here N is the total column density of the shell and n is the background number density of 
the gas to be swept up, where we assume the value for the halo gas density from model Al by 
Suchkov et al. (1994). NMgii is the column density, AMg the relative elemental abundance 
of Mg, Z/ Zq the metallicity in terms of solar abundances, and fMgii the ionization fraction 
of Mgll. The latter is quite uncertain. Z/Zq is probably between 0.1 and 1, judging from 
the abundances of damped Lya systems at low z (e.g. Boisse et al 1998). The radius R 
is again fully consistent with observed superbubbles (if much of the interstellar medium is 
swept up as well the radius could be smaller, of course). 

It is conceivable though that the z — 0.5656 system is not produced in the bubble walls 
of a smallish (~ 1 kpc) shell but rather by superwinds. Then the absorption would come 
from the dense shreds and filaments formed when the interstellar medium gets entrained by 
the hot outflowing wind (cf Bond et al 2001; Suchkov et al 1996; Heckman et al 2000, 2001; 
Mori, Ferrara & Madau 2002). The ragged component structure seen here may be caused 
by differential motion between gas clouds accelerated by the surrounding wind. 

The possibility of the absorber being caused by (symmetric ?) infall of gas has also 
been discussed, e.g., in the form of high velocity clouds which may be sweeping up galactic 
material. Apparently sizes and infall velocities of known high velocity clouds are sufficient 
to produce large scale features in the galactic gas upon impact and may even be able to 
more easily satisfy the energetic requirements of the largest shells than collective supernova 
explosions (e.g., Oort 1967; Tenorio-Tagle 1981; Tenorio-Tagle & Bodenheimer 1988 and 
references therein; Rand & Stone 1996). At high redshift low ionization metal absorption 
systems from infalling gas are expected to arise as a natural by-product of galaxy formation 
(Rauch, Haehnelt & Steinmetz 1997). Once the gas makes it past the accretion shock 
surrounding the forming galaxies it cools and can be seen in Mgll or damped Lya absorption. 
However, it is not clear whether cosmological infall of pre-enriched material can reproduce all 
the small scale density and velocity structure seen in real metal absorption systems (paper 
I, II). 

On the basis of the observations presented here and the rather idealized theoretical 
models we can currently not elucidate the origin of the z=0.5656 system any further. It would 
be interesting to study more Mgll systems in the spectra of lensed QSOs to see whether there 
is a statistical difference in spatial and velocity extent between the two-component systems 
like the one seen here and randomly selected systems. 
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6. The z=0.82743 system 

This system also has a prominent double-component structure that appears in all lines of 
sight, albeit shifted and with rather different column densities (Figs. 7 and 8). The velocity 
separation Av between the two most prominent Mgll lines in each of the A, C and B images 
is 35.0, 47.9 and 18.6 kms~^, respectively. The velocity structure is more complicated than 
in the z=0.5656 system in that here are also shifts of the column-density weighted velocity 
along the line of sight, v, between the lines of sight, which amount to vc — vX — —32.2 
kms"-*^, vb — va — 16.5 kms"-*^, and vb — vc — 64.2 kms"-*^. 

7. The z=0.97163 system 

For this system (Figs. 9 & 10) only Fell hnes are covered. Again there is some resem- 
blance between the absorption patterns for the three hnes of sight. The differences between 
the column density weighted velocities along the line of sight are vc — va — 4.1 kms"^^, 
vb—va = —13.4 kms""^, and vb — vc = —17.4 kms~^. The differential motion in this and the 
z=0.8274 system are small enough to be attributed to rotation, .e.g., of an entire gaseous halo 
(e.g., Lanzetta & Bowen 1992; Steidel et al 2002), but even for the simplest rotation model 
there are too many free parameters to be constrained by the anecdotal evidence presented 
here, and a larger sample is clearly required. 

8. Conclusions 

(1) Call absorption from the Galaxy is found in all three lines of sight to Q2237+0305 
with an average heliocentric velocity —13.9 ± 0.5 kms^^ and without significant differences 
between the lines of sight. This velocity is consistent with the absorption arising in the local 
bulk flow in the ISM within ~ 30 pc from the sun (Frisch et al 2002). In that case the largest 
angular separation between the sighthnes corresponds to < 1.8" x 30 pc = 54 astronomical 
units. Apparently, the Galactic Call cloud seen toward Q2237-I-0305 is smooth on scales 
smaller than that. 

(2) We have detected a high velocity cloud system related to the lensing galaxy 2237-1-0305 
(z=0.039). Unlike Call systems in our Galaxy this one does not show a strong component at 
the systemic velocity of the galaxy as measured either from the stellar or HI centroid. The 
system arises within impact parameters on the order of 1 kpc to the center of the galaxy (i.e., 
behind and / or in front of its bulge), so the large velocity differences among the components 
along individual lines of sight and between different lines of sight cannot be a result of the 
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rotation of the galactic disk. The absorption hne strengths are similar to those measured for 
intermediate and high velocity clouds in our Galaxy. The system may have an origin in the 
bulge or halo of the galaxy 2237+0305. 

(3) Three Mgll or Fell selected low ionization systems not associated with the lensing 
galaxy show coherent absorption over the transverse separation between the lines of sight, 
ranging from about 200 - 600 h^^^ pc. This is a strict lower limit to the size of the the gaseous 
structures giving rise to the entire Mgll system: they must be at least that large to cover all 
three lines of sight. The closeness of the lines of sight in the present case does not permit 
an upper limit to be obtained. 

(4) individual Mgll components are much smaller than entire Mgll systems: interpreting 
the separation where the optical depth between two lines of sight differ by 50% as a minimum 
'size' of the clouds, the three systems examined here show fluctuations at that level for 
separations at least as small as 200 — SOOh^Q pc. In fact, several components do not show 
up at all in more than one image (Figs. 5-10). This is consistent with the Petitjean et al 
upper limit but it indicates that many of the cloudlets are actually much smaller than a kpc. 
This result appears in contrast to size estimates for CIV systems. Individual CIV clouds 
appear to be larger than about SOOh^Q pc (paper II). Ionization calculations place the Mgll 
gas at a higher density than the more common, higher ionization, CIV gas clouds. Thus it 
is plausibly that the Mgll cloudlets are smaller than those clouds as well. 

(5) We have tentatively interpreted the low ionization system at z — 0.5656 as a moving 
shell, possibly caused by an expanding super bubble or galactic wind. The large scale coher- 
ence of the system and the differential velocities between the absorption in the three lines of 
sight are consistent with supershell parameters similar to those observed in the Galaxy and 
in local dwarf galaxies. In this picture the two Mgll components observed may correspond 
to the wall of the cooling shell at the interface with the intergalactic medium and / or the 
entrained material within the hot bubble. Observations of transitions associated with higher 
ionization stages, especially CIV and OVI, in the UV should be able to shed further light on 
the nature of this system. 

We thank Crystal Martin for helpful discussions, and Bob Carswell and Gary Ferland 
for making VPFIT and CLOUDY, respectively, available. At Keck, Gary Puniwai, Ron 
Quick, Gabrelle Saurage, David Sprayberry, Wayne Wack, and Greg Wirth helped with 
the observations. The anonymous referee is thanked for several helpful suggestions. MR is 
grateful to the NSF for grant AST-0098492 and to NASA for grant AR-90213.01-A. The 
work of WLWS was supported by NSF grant AST-9900733. 
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9. Appendix 

This is an outline of the calculation of the relation between the expansion velocity v of 
a spherical shell, and its radius i?, as a function of the observed velocity spreads between 
absorption line pairs like the one in the z = 0.5656 Mgll system {v^ = 69.8 kms~^, vb = 102.5 
kms~^, vc — 81.4 kms"-*^) and the geometric constraints (the projected relative positions of 
the images A, B, and C, as given by the distances AB (=0.66 kpc) and BC (=0.50 kpc), 
and the angle Z(ASC)(= 47.5 deg)). 

There are four unknowns, which can be expressed as (for example) i?, h (the impact 
parameter of line of sight B with respect to the center of the sphere), and (i.e., the 
angle between h and AB, the line connecting lines of sight A and B), and three equations 
connecting them (see below). Thus we can only solve for a relation between any two of these. 
The relation v{R) is probably the most useful one as it determines which values of size and 
expansion velocity are consisting with the observations. The various geometric relations 
between these quantities can be written as: 



«Ml-(|f)>.^ (3) 



(^1 - j =b' + AB^ + 2hAB cos (4) 

{l - (^)') = b^ + BC^ + 2bBCcos{(f) + Z{ABC)) (5) 
Elimination of the unwanted parameters leads to the single equation 



R-^-SAB-" BC v^A BC + R\vl - 4)) y+ 
2AB BCR'ivl - vl){4BC\^ + R^{vl - vl))y+ 
16AB^BC\^ + BC^R\vl - vlf + Ab\iQB&v^ + R\vI - vlf- 

8BC^R%\-v^c-^W + 8vh^-v\z^ + vlil + y^-z^)))) = 0, (6) 

where y = cos{Z{ABC)), and z = sm{Z{ABC)). It can be solved for the velocity v. Of 
the four solutions, two are imaginary and one describes a contracting {v < 0) rather than 
an expanding sphere. The remaining solution is well-behaved for 0.33 kpc < R < oo: 



V ^ 1/2 [(AB BC R\v\ - v%)y + AB BCR\v'c - v%)y + 
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[{Ab'bc'r\-{BC'' - 2AB BCy - ARh'' + Ab')(Ab'(4 - v^f - 
2AB BC{vl - vDivl - vl)y + BC\vl - vlf) + 

(BC\vl-v^y + A^ivl- vl)y + AB^{vl + vl{z' - 1 - y^) + v\)f)Y'')l 
(AB^BC'^(BC^ - 2AB BCy - AR'^z'^ + AB^))]^/^ (7) 



Inserting the numerical values gives 



0.5 



4.878 X 1015. 



i?4(0.3959 + 



+ 



3.4712 X 10i5i?2 



(9.3456 - 84.5541i?2)2 (84.5541i?2 - 9.3456) ' 



(8) 



where R is now to be measured in units of kpc, and v in cms ^. The relation v{R) is shown 
in Fig. 11. 



-19- 



REFERENCES 

Bergeron, J., Boisse, P., 1991, A&A, 243, 344 

Bergeron, J.; Stasinska, G., 1986, A&A, 169, 1 

Blades, J.C., Morton, D.C., 1983, MNRAS, 204, 317 

Boisse, P., Le Brun, V., Bergeron, J., Deharveng, J.-M., 1998, a, 333, 841 

Bond, N.A., Churchill, C.W., Charlton, J.C., Vogt, S.S. 2001, ApJ, 557,761 

Bond, N.A., Churchill, C.W., Charlton, J.C., Vogt, S.S. 2001, ApJ, 562,641 

Bowen, D.V., 1991, MNRAS, 251, 649 

Bowen, D.V., Roth, K.C., Blades, J.C., Meyer, D.M., 1994, ApJ,420, 71 
Bowen, D.V., Blades, J.C., Pettini, M., 1995, ApJ,448, 634 

Carswell, R.F.,Webb, J. K., Cooke, A.J., Irwin, M.J., 1992, Institute of Astronomy, Cam- 
bridge, http://www.ast.cam.ac.uk/ rfc/vpfit.html 

Carswell, R.F.,Webb, J.K., Lanzetta, K.M., Parnell, H.C., Webb, J.K., 1991, ApJ, 371, 36 

Churchill, C.W., Rigby, J.R., Charlton, J.C., Vogt, S.S., 1999, ApJS, 120, 51 

Crane, P., Albrecht, R., Barbieri, C, Blades, J.C., Bokensenberg, A., Deharveng, J.M., 
Disney, M.J., Jakobsen, P., Kamperman, T.M., King, I.R., Macchetto, F., Mackay, 
CD., Paresce, F., Weigelt, G., Baxter, D., Greenfield, P., Jedrzejewski, R., Nota, A., 
Sparks, W.B., 1991, ApJ, 369, 59 

Deul, E.R., den Hartog, R.H., 1990, A&A, 229, 362 

D'Odorico, S., di Serego Alighieri, S., Pettini, M., Magain, P., Nissen, P. E., Panagia, N., 
1989, A&A, 215, 21 

Ferland, G.J., 1993, University of Kentucky Department of Physics and Astronomy Internal 
Report 

Foltz, C.B., Hewett, P.C., Webster, R.L., Lewis, G.F., 1992, ApJ, 386, L43 
Frisch, P.C., Grodnicki, L., Welty, D.E., 2002, ApJ, submitted (astro-ph/0203093) 
Greenstein, J.L., 968, ApJ,152, 431 

Heckman, T.M., Lehnert, M.D., Strickland, D.K., Armus, L., 2000, ApJS, 129, 493 

Heckman, T.M., Sembach, K.R., Meurer, G.R., Strickland, D.K., Martin, C.L., Calzetti, D., 
Leitherer, C., 2001, ApJ, 554, 1021 

Heiles, C, 1979, ApJ, 229, 533 

Heiles, C, 1976, ApJ, 207, L137 



-20- 



Hintzen, P., Maran, S.P., Michalitsianos, A.G., Foltz, C.B., Chaffee, F.H., Kafatos, M., 1990, 
AJ, 99, 45 

Ho, L.C., Filippenko, A.V., 1995, ApJ, 444, 165 

Hucfira, J., Gorenstein, M., Kent, S., Sfiapiro, I., Smitfi, G., Horine, E., Perley, R., 1985, 
AJ, 90,691 

Irwin, M.J., Webster, R.L., Hewett, P.C., Corrigan, R.T., Jedrzejewski, R.I., 1989, AJ, 98, 
1989 

Kennicutt, R.C., Bresolin, F., Bomans, D.J., Botliun, G.D., Tliompson, LB., 1994, AJ, 438, 
594 

Lanzetta, K.M., Bowen, D.V.,1992 ApJ, 391, 48L 

Lopez, S., Reimers, D., Rauch, M., Sargent, W.L.W., Smette, A., 1999, ApJ, 513, 598 

Meyer, D.M., Rotli, K.C., 1991, ApJ, 383, L41 

Monier, E.M., Turnsliek, D.A., Lupie, O.L., 1998 ApJ, 496, 177 

Morton, D.C., Andereck, C.D., 1976, ApJ, 205, 356 

Morton, D.C., Blades, J.C., 1986,MNRAS, 220,927 

Oey, M.S., Clarke, C.J., 1997, MNRAS, 289, 570 

Oort, J.H., 1967, lAU Symp. No. 31, 279 

Pctitjcan, P., Aracil, B., Srianand, R., Ibata, R., 2000, A&A, 359, 457 
Rand, R.J., Stone, J.M., 1996, AJ, 111 190 

Rauch, M., Sargent, W.L.W., Womble, D.S., Barlow, T.A., 1996 ApJ, 467, 5 

Rauch, M., Haehnelt, M.G., Steinmetz, M., 1997, ApJ, 481, 601 

Rauch, M., Sargent, W.L.W., Barlow, T.A., 1999, ApJ, 515, 500 (paper I) 

Rauch, M., Sargent, W.L.W., Barlow, T.A., 2001, ApJ, 554, 823 (paper II) 

Rauch, M., Sargent, W.L.W., Barlow, T.A., Carswcll, R.F., 2001, ApJ,562, 76 (paper III) 

Rigby, J.R., Charlton, J.C., Churchill, C.W., 2002, ApJ, in press. 

Robertson, J.G., Schwarz, U.J., van Woerden, H., Murray, J.D., Morton, D.C., Hulsbosch, 
A.N.M., 1991, MNRAS, 248, 508 

Smette, A., Robertson, J.G., Shaver, P.A., Reimers, D., Wisotzki, L., Koehler, T., 1995, 

A&A, 113, 199 

Songaila, A., 1981, ApJ248, 945 

Songaila, A., Blades, J.C., Hu, E.M., Cowie. L.L., 1986, ApJ, 303, 198 



- 21 - 



Tenorio-Tagle, G., 1981, A&A94, 338 

Steidel, C.C., Rich, R.M., McCarthy, J.K., 1990, A J, 99, 1476 
Steidel, C.C., 1990, ApJS, 74, 37 

Steidel, C.C., 1995, in ESO Workshop on Quasar Absorption Lines, ed. G. Meylan (Garching: 
Springer), 139 

Steidel, C.C., Kollmcicr, J.A., Shapley, A.E., Churchill, C.W., Dickinson, M., Pettini, M., 
2002, ApJ, in press. 

Suchkov, A.A., Balsara, D.S., Heckman, T.M.; Leitherer, C, 1994, ApJ, 430, 511 

Suchkov, A.A., Berman, V.G., Heckman, T.M. Balsara, D.S., 1996, ApJ, 463, 528 

Vladilo, G., Centurion, M., de Boer, K.S., King, D.L., Lipman, K., Stegert, J.S.W., Unger, 
S.W., Walton, N.A., 1995, A&A, 291, 425 

Vogt, S.S., et al., 1994, S.P.I.E. 2198, 362 

Wakker, B.R, van Woerden, H., Gibson, B.K., 1999, in ASP Conf.Ser. 166, Stromlo Work- 
shop on High Velocity Clouds, ed. B.K. Gibson & M.E. Putman, (San Francisco: 
ASP), 311 

Wakker, B.R, Mathis, J.S., 2000, ApJ, 544, 107 
Wakker, B.R, 2001, ApJS, 136, 463 

West, K.A., Pettini, M., Penston, M.V.; Blades, J.C.; Morton, D.C., 1985, MNRAS, 215, 
481 

Yee, H.K.C.. 1988. AJ. 95. 1331 



This preprint was prepared with the AAS lAT^TjX macros v5.0. 



-22- 



Tablc 1. The Galactic Call system 



LoS ion z b logN 



A Call -0.000043± 0.000002 15.21± 1.06 12.259± 0.025 
B Call -0.000056± 0.000004 15.88± 1.73 12.194± 0.039 
C Call -0.000048± 0.000003 IS.Oli 1.20 12.256± 0.024 
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Table 2. The z=0.039 (lens) system 



LoS 


ion 




h 


logA^ 


A 


Call 


0.038700± 0.000006 


8.76± 2.76 


11.513± 0.093 


A 


Call 


0.039369± 0.000002 


7.63± 1.00 


11.901± 0.039 


B 


Call 


0.038016± 0.000014 


18.78± 7.12 


11.722± 0.175 


B 


Call 


0.038316± 0.000018 


24.34± 8.91 


11.753± 0.169 


B 


Call 


0.039377± 0.000013 


12.80± 6.05 


11.490± 0.171 


C 


Call 


0.038483± 0.000006 


15.46± 2.84 


11.842± 0.072 


c 


Call 


0.038685± 0.000001 


10.18± 0.45 


12.513± 0.016 


c 


Call 


0.039146± 0.000010 


15.05± 4.63 


11.604± 0.114 


c 


Call 


0.039325± 0.000009 


17.06± 4.05 


11.753± 0.093 
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Table 3. The z=0.566 system, with corrected zero levels 



LoS 


ion 




h 

yj 


lop- N 

lUg iV 


A 


Mgll 


0.565696± 0.000003 


4.16± 1.26 


11.801± 0.059 


A 


Mgll 


0.566289± 0.000004 


8.92± 1.30 


13.020± 0.105 


A 


Mgll 


0.566304± 0.000003 


3.79± 1.74 


12.844± 0.196 


A 


Mgll 


0.566352± 0.000035 


22.91± 4.53 


12.987± 0.044 


A 


Mgll 


0.566480± 0.000038 


5.94± 5.16 


12.383± 0.829 


A 


Mgll 


0.566518± 0.000002 


4.72± 0.59 


13.484± 0.189 


A 


Mgll 


0.566655± 0.000003 


9.49± 1.16 


12.303± 0.039 


A 


Mgll 


0.566797± 0.000010 


8.32± 3.53 


11.622± 0.119 


A 


Mgll 


0.566929± 0.000001 


4.97± 0.25 


12.730± 0.022 


A 


Mgl 


0.566289± 0.000000 


8.80± 4.28 


11.329± 0.195 


A 


Mgl 


0.566304± 0.000000 


0.40± 0.78 


12.026± 2.909 


A 


Mgl 


0.566352± 0.000000 


27.50± 30.69 


10.945± 1.180 


A 


Mgl 


0.566480± 0.000000 


24.91± 17.02 


11.296± 0.425 


A 


Mgl 


0.566518± 0.000000 


5.01± 1.41 


11.459± 0.097 


A 


Mgl 


0.566655± 0.000000 


65.20± 146.68 


10.929± 1.058 


A 


Mgl 


0.566929± 0.000000 


7.23± 1.72 


11.311± 0.080 


A 


Fell 


0.565696± 0.000000 


32.19± 38.86 


11.974± 0.829 


A 


Fell 


0.566289± 0.000000 


11.14± 1.02 


13.365± 0.078 


A 


Fell 


0.566304± 0.000000 


1.72± 2.01 


12.751± 0.239 


A 


Fell 


0.566352± 0.000000 


26.25± 7.23 


12.991± 0.192 


A 


Fell 


0.566480± 0.000000 


8.54± 4.75 


12.764± 0.607 


A 


Fell 


0.566518± 0.000000 


4.76± 0.80 


13.405± 0.154 


A 


Fell 


0.566655± 0.000000 


8.30± 2.44 


12.313± 0.160 


A 


Fell 


0.566797± 0.000000 


21.01± 44.96 


11.891± 0.802 


A 


Fell 


0.566929± 0.000000 


4.80± 1.12 


12.496± 0.082 


B 


Mgll 


0.566180± 0.000002 


4.37 ± 0.58 


14.222 ± 0.472 


B 


Mgll 


0.566270± 0.000003 


8.88 ± 0.71 


13.100 ± 0.034 


B 


Mgll 


0.566432± 0.000003 


4.67 ± 0.85 


12.311 ± 0.057 


B 


Mgll 


0.566496± 0.000001 


3.81 ± 0.80 


12.627 ± 0.037 


B 


Mgll 


0.566570± 0.000002 


4.02 ± 0.93 


12.469 ± 0.039 


B 


Mgll 


0.566660± 0.000003 


6.15 ± 0.86 


12.852 ± 0.055 


B 


Mgll 


0.566709± 0.000015 


5.27 ± 4.07 


11.776 ± 0.455 


B 


Mgll 


0.566825± 0.000007 


7.79 ± 2.90 


11.869 ± 0.101 
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Table 3 — Continued 
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Table 3 — Continued 
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z 


h 


log TV 


C 


Fell 


0.566607± 0.000000 


4.76± 0.90 


13.085± 0.045 


c 


Fell 


0.566871± 0.000000 


0.18± 0.32 


12.019± 0.727 
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Table 4. The z=0.827 system 
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z 


h 


logiV 


A 


Mgll 


0.827437 ± 0.000001 


2.53±0.18 


12.585± 0.039 


A 


Mgll 


0.827638 ± 0.000004 


1.63±2.17 


11.520± 0.214 


A 


Mgll 


0.827682 ± 0.000060 


7.54±12.58 


11.112± 0.644 
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Fell 


0.827435 ± 0.000002 
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B 


Mgll 


0.827646± 0.000004 
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B 


Mgll 
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1.33±0.24 
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B 
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B 


Fell 


0.827795± 0.000003 
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12.089± 0.146 


C 


Mgll 


0.827461± 0.000008 


2.90 ±3.27 
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C 
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0.827608± 0.000004 


9.59 ±1.14 


11.923± 0.037 
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Table 5. The z=0.972 system 



LoS 


ion 


z 


h 


logA^ 


A 


Fell 


0.971635± 0.000001 


3.35± 0.23 


13.362± 0.021 


A 


Fell 


0.971668± 0.000010 


8.99± 0.94 


12.570± 0.119 


A 


Fell 


0.971965± 0.000018 


4.06± 6.21 


11.266± 0.265 


A 


Fell 


0.971866± 0.000020 


3.42± 8.19 


11.113± 0.361 


B 


Fell 


0.971545± 0.000002 


4.02 ± 0.47 


12.505± 0.025 


B 


Fell 


0.971646± 0.000019 


6.57 ± 5.81 


11.579± 0.211 


C 


Fell 


0.971576± 0.000008 


4.95 ±1.49 


12.153 ± 0.246 


C 


Fell 


0.971651± 0.000041 


9.44 ±6.83 


12.028 ± 0.333 


C 


Fell 


0.972003± 0.000007 


6.17 ±1.96 


11.778 ± 0.084 




Fig. 1. — Keck HIRES guider camera image of Q2237+0305. The exposure time was 5 s, the 
picture was taken on October 21, 1998. The relative flux ratios at that date were approximately 
Fa : Fb : Fc : Fd = 1.0 : 2.1 : 1.7 : 0.7. 
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Call 3970 (B vs. C) 



Call 3935 (B vs. C) 



Call 3970 (A vs. C) 



Call 3935 (A vs. C) 




Call 3970 (A vs. B) 



Call 3935 (A vs. B) 



-200 200 

heliocentric velocity [kms ''] 

Fig. 2. — absorption spectra of the Call 3935, 3970 A absorption from the Milky Way (hehocentric 
velocities). The three lines of sight A, B, and C, arc compared in pairs. The uppermost three spectra 
show the two transitions in lines of sight B (thin line) versus C (thick line), the middle two spectra 
show A (thin line) versus C (thick line), and the bottom two spectra A (thin line) versus B (thick 
line), respectively. Here and in the other figures the spectra are normalized to a unit continuum 
and offset along the ordinate (labelled in abitrary units) for better visibility. 
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Call 3970 (B vs. C) 



Call 3935 (B vs. C) 



Call 3970 (A vs. C) 




Call 3935 (A vs. C) 



Call 3970 (A vs. B) 



Call 3935 (A vs. B) 



-200 



200 



velocity [kms ^] relative to z = 0.039014 

Fig. 3. — spectra of the Call 3935, 3970 A arising in the lensing galaxy (the origin of the velocity 
scale corresponds to the galaxy rcdshift given by Foltz et al. 1992). The three lines of sight A, B, 
and C, are again compared in pairs. The lines marked with "x" are interlopers of metal lines at 
different redshift. The proper beam separation between A and B is here 1.89 h^Q kpc, between A 
and C 1.41 h7} kpc, and between B and C 1.45 hz} kpc. 
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-200 200 

velocity [kms ] relative to z = 0.039014 

Fig. 4. — Same system as in Fig. 3, (Call at z=0.038) but now showing only the 3935 A transition 
with the three systems overplotted offset by the height of the normahzed continuum. The three 
hues of sight A, C, and B, are arranged from bottom to top such that the separation from image A 
increases as one goes from the bottom to the top spectrum. The proper beam separation between 
A and B is here 1.89 h'^^ kpc, between A and C 1.41 h^^ kpc, and between B and C 1.45 h'^^ kpc. 
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200 400 

velocity [kms ^] relative to z = 0.56561 1 



Fig. 5. — absorption spectra of tlie tliree transitions Mgl 2853, Fell 2600, and Mgll 2796 in the 
strong low ionization system near z=0.5656. The three hnes of sight A, B, and C, are compared in 
pairs. The uppermost three spectra show the three transitions in hnes of sight B (thin hne) versus 
C (thick hne), the middle three spectra show A (thin line) versus C (thick line), and the bottom 
three spectra A (thin line) versus B (thick line), respectively. The proper beam separation between 
A and B is here 0.66 h^Q kpc, between A and C 0.49 h'^^ kpc, and between B and C 0.50 h^Q kpc. 
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100 200 

velocity [kms ] relative to z = 0.56561 1 

Fig. 6. — absorption spectra of the transition Mgll 2796 (at z=0. 56561), in the strong low ionization 
system shown already in the previous figure. The three lines of sight A, B, and C are overplotted 
with offsets of 1.0 for clarity. There is a minimum of absorption near 150 kms~^ on the arbitrary 
velocity scale. Also, going from the bottom to the top spectrum the separation from the A image 
increases (moving in a north-east direction on the sky - see Fig. 1) and the two stronger groups 
separated by the 150 kms^^ gap move further apart. The proper beam separation between A and 
B is 0.66 h^Q kpc, between A and C 0.49 h^Q kpc, and between B and C 0.50 h^Q kpc. 
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Fell 2600 (B vs. C) 



Fell 2600 (A vs. C) 
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200 400 

velocity [kms ^] relative to z = 0.82743 

Fig. 7. — absorption spectra of the two transitions Fell 2600, and Mgll 2796 in the strong low 
ionization system near z=0.827. The three lines of sight A, B, and C, are compared in pairs. The 
uppermost two spectra show the two transitions in lines of sight B (thin line) versus C (thick line), 
the middle two spectra show A (thin line) versus C (thick line), and the bottom three spectra A 
(thin line) versus B (thick line), respectively. The proper beam separation between A and B is 0.38 
h'^Q kpc, between A and C 0.29 h'^Q kpc, and between B and C 0.30 h'^^ kpc. 
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Fig. 8. — Same system as in Fig. 7 (z=0.827), but now showing only Mgll 2796. All three systems 
are overplottcd offset by the height of the normalized continuum. The three lines of sight A, C, 
and B, are arranged from bottom to top such that the separation from image A increases as one 
goes from the bottom to the top spectrum. The proper beam separation between A and B is 0.38 
h^Q kpc, between A and C 0.29 h'^Q kpc, and between B and C 0.30 h^Q kpc. 
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velocity [kms ] relative to z = 0.97163 

Fig. 9. — absorption spectra of the Fell 2383 transition low ionization system near z=0.97163. 
The three Hnes of sight A, B, and C, are compared in pairs. The uppermost spectrum show the 
transition in lines of sight B (thin line) versus C (thick line), the middle spectrum show A (thin 
line) versus C (thick line), and the bottom spectrum A (thin line) versus B (thick line), respectively. 
The proper beam separation between A and B is 0.28 h'^Q kpc, between A and C 0.21 h'^^ kpc, 
and between B and C 0.22 h^Q kpc. 
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Fig. 10. — Same system as in Fig. 9, (Fell 2353 at z=0. 97163) but now with the all three systems 
overplotted offset by the height of the normalized continuum. The three lines of sight A, C, and 
B, are arranged from bottom to top such that the separation from image A increases as one goes 
from the bottom to the top spectrum. The proper beam separation between A and B is 0.28 ^^q^ 
kpc, between A and C 0.21 h^^ kpc, and between B and C 0.22 h^Q kpc. 
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Fig. 11. — expansion velocity of a spherical bubble consistent with our observations, as a function 
of radius. The curve v{R) was computed from equation 8. For comparison the hatched area shows 
the approximate ranges of sizes and expansion velocities of HI shells and supershells from the 
literature. 
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Fig. 12.— 



three lines of sight to QSO 

schematic drawing of an expanding spherical bubble intersected by three lines of sight. 



